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I. INTRODUCTION 


The primary objective of the present follow on to Contract 
NAS8-26131 is the addition of a dynamic loads computation capability to 
the Unified Flexible Spacecraft Simulation (UFSS) Program originally 
developed for the NASA MSFC under the above named contract. This added 
capability provides a means for determining the internal member loads due 
to the time-varying external loading conditions experienced by an orbiting 
spacecraft/cluster. 

Consistent with the guidelines used in the development of the 
original UFSS Program, the dynamic loads computation is completely 
generalized and applicable to structural components of widely varying 
configurations . 

As described below, the present updated UFSS Program and the 
newly developed Unified Flexible Spacecraft Load Program (LOAD) completely 
fulfill all the objectives of the original and follow-on contracts. 
Following the initial delivery to the MSFC, several efficiency modifica- 
tions and advanced capabilities were added to the UFSS Program by TRW 
Systems in the normal course of using the program in-house. These 
capabilities (to be described in Section IV below) have been included in 
the present version of the program at no extra cost to the MSFC. 

A complete list of final deliverables for the present contract is 
contained in Section III of this document. 

II. GENERAL DESCRIPTION OF THE DYNAMIC LOAD, CALCULATION CAPABILITY 

The UFSS Program is a general purpose digital computer program 
designed to simulate the dynamic response of controlled, multibodied, 
flexible space vehicles subjected to environmental and operational 
disturbances. The synthesizing algorithm within. UFSSP is based on elimi- 
nation of the interaction forces and torques between adjacent bodies of 
the system model. By so doing, the final system of matrix differential 
equations includes only the unconstrained degrees of freedom of the 
spacecraft model. This retention of the minimal number of degrees of 
freedom contributes greatly to reducing the cost of the solution for the 

However, when the dynamic loads are desired, the 


dynamic response. 



interaction forces and torques, as well as the internal loads within a 
given flexible body, must be obtained explicitly. 

Determination of these loads can be divided into two basic phases. 

In the first phase, the forces and moments acting at each interconnection 
between adjacent bodies are calculated within the UFSS Program via 
information available from the original dynamics subroutines. In the 
second phase, the mode-acceleration method is utilized to calculate 
internal loads within any desired member of the structural model for a 
given terminal flexible body. These internal loads are calculated by a 
separate, stand-alone program operating upon a special loads history tape 
generated by the UFSS Program. Thus, if only the interconnection loads 
are desired, these can be obtained from the UFSS Program time history run 
itself by merely setting the LOAD option flag to "ON". User information 
necessary to run the UFSS Program is contained in Reference 1, while 
details of the equations actually programmed to perform these load calcu- 
lations are presented in Appendix D of Reference 2. Programming operation 
of the UFSS Program is contained in Reference 3. 

Because of the core storage required to save the necessary input 
data and the fact that internal load calculations can be performed sub- 
sequent to the dynamic simulation of a given system, it was decided that 
the stand-alone program "LOAD" be generated to calculate the internal 
loads within a given flexible body instead of calculating them within a 
new subroutine of UFSSP. 

In the UFSS Program, flexible bodies are modeled in the traditional 
structural dynamics sense as a system of joints (or nodes) which are 
interconnected by weightless finite element members (e.g., beams, plates...). 
All masses are lumped at the joints. The orthogonal functions used to 
describe the spatial deformation of the bodies are normally taken to be the 
orthonormal cantilever modes produced by a standard structural dynamics 
program such as NASTRAN, SAMIS, or SMAP. In general, most such programs 
are based on small deflection theory, using the direct stiffness matrix 
finite element approach assuming linear stiffnesses. Such an approach 
allows for the generation of a specific transformation matrix, herein 
called the load transformation matrix (LTM) , by a systematic application of 
a unit force along each degree of freedom with all other forces set equal 
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to zero. In particular, such a matrix allows the LOAD Program to calculate 
internal member forces through the simple matrix multiplication operation: 


V c) ■ b mq F q (t) 


( 1 ) 


In the above equation, each component of L^ (t) represents a specific 
desired internal member load, is the LTM and ( t ) is the applied force 

vector. Specifically, contains, as rows, coefficients for each degree 

of freedom relating the desired load to a unit force applied successively 
along each translational degree of freedom of the structure (assumed to 
be three times the number of joints in the UFSSP model since only the modal 
translations are used to describe deformations), and the ordering of the 
columns must be identical to the degree of freedom ordering within the 
mode shapes. The applied force vector, F^(t), contains both the externally 
applied forces (environmental forces, control forces...) and the inertial 
(d' Alembert) forces associated with the inertial velocity and acceleration 
of each joint. For the present initial implementation of the LOAD Program, 
the allowable externally- applied forces acting on the flexible bodies 
can be added at any future time if a specific application calls for their 
inclusion), there are no restrictions on the disturbance forces and/or 
torques acting on the rigid bodies of the system model. 

Figure 1 presents an overview on the program interfaces. As seen 
from Equation (1), only two basic quantities are required for the loads 
computation. The first quantity, B^ , is obtained from any suitable 
structural dynamics program; the second quantity, F (t) , is obtained 
from the UFSS Program via a special history ' tape . 

1 / 
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Figure 1. Basic Interfaces for Internal Load Calculations 


User information necessary to run the LOAD Program is contained 
in Reference 4, while details of the equations actually programmed to 
perform the LOAD calculations are presented in Appendix D of Reference 2. 
Since the LTM must be obtained from an independent program of the user's 
choice, descriptions of its formulation and format are contained in 
Section III and Appendix A of Reference 4. 
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III. LIST OF DELIVERABLES UNDER THE PRESENT CONTRACT 


The following items constitute the final deliverables for the 
present follow-on to Contract NAS 8-26131: 

1. A PCF tape containing the source and object decks for the 
updated UFSS Program suitable for operation on the MSFC 
Univac 1108 EXEC II computing system. 

2. A PCF tape containing the source and object decks for the 
LOAD Program suitable for operation on the MSFC Univac 1108 
EXEC II computing system. 

3. Input data decks necessary to run the UFSS Program "user 
manual test case" plus print-out of the test case results. 

4. Input data decks and a copy of the UFSSP history tape necessary 
to run the LOAD Program "user and operation manual test case" 
plus print-out of the' test case results. 

5. Nine (9) copies of this report and each of the documents listed 
in the References section of this report. 

IV. ADDITIONAL PROGRAM FEATURES INCLUDED IN THE UFSS PROGRAM 

In the course of using the UFSS Program in-house, the following 
features have been incorporated into the UFSS Program and are hereby 
included gratis in the present delivery of the program to the MSFC: 

1. In addition to the original capability of inputting the flexible 
structures data via magnetic tape, a simplified format has 

been included for inputting this data via a card deck or a 
disk/drum file. 

2. A "high spin" capability has been included in the program. 

The dynamics of rapidly spinning flexible bodies (modeled 

as either space curves or flat plates) is now handled automatically 
by an internal addition of a modified displacement function 
(curvature shortening effect) for each flexible body; thus, 
the centrifugal stiffening effect as well as all other 
acceleration effects are now included in the flexible body 
dynamics. \J ( t 
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In the initial UFSS Program, Body 1 kinematics were obtained 
by direct integration of its nine element direction cosine 
matrix. In the present version, this method has been replaced 
by a much more efficient and accurate calculation of the 
direction cosine matrix via integration of four "Euler para- 
meters". 
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INTRODUCTION 


The LOAD Program is a digital computer program designed to calculate 
internal loads within the flexible bodies of any model admissible to the 
Unified Flexible Spacecraft Simulation (UFSS) Program. 

A 

As described in [1] , the UFSS Program is a general purpose digital 
computer program designed to simulate the dynamic response of controlled, 
multibodied, flexible space vehicles subjected to environmental and operational 
disturbances. The synthesizing algorithm within UFSSP is based on elimination 
of the interaction forces and torques between adjacent bodies of the system 
model. By so doing, the final system of matrix differential equations in- 
cludes only the unconstrained degrees of freedom of the spacecraft model. 

This retention of the minimal number of degrees of freedom contributes 
greatly to reducing the cost of the solution for the dynamic response. 

However, when the dynamic loads are desired, the interaction forces and 
torques, as well as the internal loads within a given flexible body, must 
be obtained explicitly. 

The forces and moments acting at each interconnection between 
adjacent bodies are calculated within the UFSS Program itself via information 
available from the original dynamics subroutines. In the LOAD program 
described herein, the mode-acceleration method is utilized to calculate internal 
loads at any desired node point within a given terminal flexible body. These 
internal loads are calculated by suitably operating upon a special loads 
history tape generated by the UFSS Program. 

This present document contains all the necessary information for 
using the LOAD program. For details of the equations utilized in the 
programming, ^the user is referred to Appendix D of the UFSSP Methodology 
Report [2], 

Section I herein contains a general description of the program and 
Section II lists the principal capabilities. Section III presents a 
complete listing and description of all input data quantities while Section 
IV contains additional input information. Section V details the available 
output from the program and the use of the error messages. Section VI 
presents sample case runs including output prints and plots. 

A 

numbers in brackets denote references. 
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Sections VII through XI constitute the Operation Manual for the 
LOAD program. Section VII contains the software/hardware requirements 
and Section VIII presents the overall deck setup. Section IX details 
the tape/disk/drum usage while Section X contains load maps for both the 
CDC 6500 and the Univac 1108. Finally, Section XI presents an overall 
block diagram of the LOAD program, while Appendix A contains a description 
of the required load transformation matrix (LTM) . 



I GENERAL DESCRIPTION OF THE LOAD PROGRAM 

In the UFSS Program, flexible bodies are modeled in the traditional 
structural dynamics sense as a system of joints (or nodes) which are 
interconnected by weightless finite element members (e.g., beams, plates...). 
All masses are lumped at the joints. The orthogonal functions used to 
describe the spatial deformation of the bodies are normally taken to be the 
orthonormal cantilever modes produced by a standard structural dynamics 
program such as NASTRAN, SAMIS, or SMAP. In general, most such programs 
are based on small deflection theory, using the direct stiffness matrix 
finite element approach assuming linear stiffnesses. Such an approach 
allows for the generation of a specific transformation matrix, herein 
called the load transformation matrix (LTM) , by a systematic application of 
a unit force along each degree of freedom with all other forces set equal 
to zero. In particular, such a matrix allows internal member forces to 
be calculated through the simple matrix multiplication operation. 

I^Ct) = Fq(t) 

In the above equation, each component of 1^(0 represents a specific 
desired internal member load, B^ is the LTM and Fq(t) is the applied force 
vector. Specifically, contains, as rows, coefficients for each degree 

of freedom relating the desired load to a unit force applied successively 
along each translational degree of freedom of the structure (assumed to be 
three times the number of joints in the UFSSP model since only the modal 
translations are used to describe deformations) , and the ordering of the 
columns must be identical to the degree of freedom ordering within the 
mode shapes. The applied force vector, Fg(t), contains both the externally 
applied forces (environmental forces, control forces...) and the inertial 
(d'Alembert) forces associated with the inertial velocity and acceleration 
of each joint. 

Because of the core storage required to save the LTM and the fact 
that internal load calculations can be performed subsequent to the dynamic 
simulation of a given system, it was decided that the present stand-alone 
program "L0AD" be generated to calculate the internal loads within a given, 
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flexible body instead of calculating them within a new subroutine of UFSSP. 
Figure 1 presents an overview of the program interfaces. As seen from 
Equation (1), only two basic quantities are required for the loads computation. 
The first quantity, B , is obtained from a structural dynamics program; 
the second quantity, F^Ct) » is obtained from the UFSS Program via a special 
history tape. 

As previously mentioned, the row index, M , of runs one 

to A, where A is the total number, of internal loads to be calculated; the 
column index, Q , runs from one to 3N, where N is the total number of 
joints in the structural dynamics model of the given flexible body, as input 
to UFSSP. Additional description of the LTM is contained in the description 
of the Sample Case presented in Section VI of this document and in Appendix A. 



Figure 1. Basic Interfaces for Internal Load Calculations 
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The LOAD Program calculates internal loads for one flexible body 
at a time; thus, a single LTM is input for each run. The input LTM is matched 
to the appropriate flexible body through the input quantity N0B0D, and the 
LTM data can be input either as a card deck or as a tape (or disk) file. 

The input LTM can contain up to 72 rows (each row corresponding to 
a desired internal load), and any subset of these rows can be selected for 
computation on a given run through input specification via the MSELEC array. 

Loads can be calculated for every time point on the UFSS history tape 
or only for specific time intervals. Through the input array STIM__n__, the 
user specifies both the time intervals for calculation and the "sampling 
rate" at which the calculations are performed. At each time point for 
which loads are calculated, all loads are written on a history tape for 
later printing and/or plotting (no printing or plotting option is incorpor- 
ated within the LOAD program itself) . 

As the values of the maximum and minimum for each load are of 
primary interest, the' standard output of the LOAD Program consists of an 
output print of these maximum and minimum values for each load, the time 
at which they occurred, and the values of all the loads at each of these 
critical times. 
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PROGRAM CAPABILITIES 


A. This program will process one case per job. This means that the 
loads will be computed for one body at a time. 

B. Plotting of Data on the history tape cannot be done in the same 
job. The TRW utility program PL0T8 can be used to plot the data 
at a later time. 

C. There can be a maximum of 4000 data time points on the output 
history tape. 

D. There can be a maximum of 72 rows of the B^. rt matrix read in core. 

MO 



III. 


LIST OF INPUT DATA ON CARDS 


The Namelist data input format is used for all card input except for 
the header. See Section IVB for instructions on how to input certain 
data items. 


A. DEFINITION OF INPUT DATA ON CARDS 


Program 

Input 

Name 

Dimension 

Nominal 

Value 

Data Description 

N0CASE 

1 

— 

8 digit case no. to identify the results 
of this job. 

N0B0D 

1 

0 

Body no. j identifying the body within 
the model processed by UFSSP. 

INPTRN 

1 

25 

Unit no. for the input of the card 
data for the loads transformation 
matrix, =5 - card read; =25 - file 
on disk or tape. 

MSELEC 

2 x 50 

All M rows 
will be 
processed 

r0W se ^ ector tables. There is a 
possibility of 50 from-to entry 
pairs. The column represents pair no. 

Row 1 is the "from" row no. and row 2 
is the "to" row no. Example: MSELEC (1,4) 
=7, 10 means select 7th through 10th 
rows of matrix from loads trans- 
formation file as the 4^ from-to 

entry pair. 

HISCSE 

1 

— 

Gives the case no. of the UFSSP history 
tape. 

STIM01- 

STIM10 

3 

0 

A set of up to 10 different time 
periods and data selection frequencies 
used to select time history data from 
the UFSSP history tape. Example: 

STIM03 (1)= . 05 , 1.0, 2 means starting 
at or after time =.05 sec. and ending 
at or before time =1 sec., select data 
for every other time point from the tape. 

INPRNT 

1 

9^0 

Print input card data option flag, 
= 0 - no; ^0 - yes. 
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IV. 


INPUT DATA SET DESCRIPTION 


A. SOURCES OF DATA OTHER THAN NAMELIST INPUT 

1. UFSSP History Tape 

This tape is generated by the UFSSP program. It contains struc- 
tures data and the necessary time history data to compute, 
internal loads. See the UFSSP User’s Manual Appendices A-D 
for the file format and Section IVC.5. for contents of a data 
record. 

2. Loads Transformation Matrix Data (B^) 

This input data is submitted either via the card reader or a 
BCD card image file on disk or tape. This data is generated 
by the LTMP program (TRW) or keypunched from data load sheets. 

Format : 

The first card contains the number of rows and columns of the B m q 
matrix which is on the cards that follow. The last data card 
is followed by an end of file mark. All cards contain a i.our 
digit sequence number in c.c. 77-80 starting with 0001 and 
increasing by -l's. 

Card 1 5 Header 

C.C. 1-3 "M" - No. of rows in matrix 

C.C. 4-6 "Q" - No. of columns in B^ 0 matrix 3 3*NJTS. 

Card 2 - Card n: Matrix 

There are up to 5 elements per card 

C.C. 1-13 Field 1; C.C. 14-26 Field 2 
C.C. 27-39 Field 3; C.C. 40-52 Field 4 
C.C. 53-65 Field, 5 

Each data field will be of the form: +X . XXXXXXE4YY 

The matrix is on the cards by rows so that all columns for one 
row will come in before the next row is read. Each row starts 
on a new card. 
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B. RULES FOR INPUTTING DATA ON CARDS 

1. General 

a. Floating Point Numbers 

+XXXXX.X, +X.XXXXXXE+YY and XXXX are all valid ways of 
inputting floating point numbers. A maximum of 8 signifi- 
cant digits is allowed. 

b . Fixed Point Numbers 

+XXXXX; up to 10 digits are allowed with no decimal point. 

c. Omitted Data 

Except for data with nominal values and the header which 
is stored as blanks, the missing items will be set to zero. 

d. Namelist Rules 

The rules for entering data under the Namelist input routine 
are .in effect. See CDC FORTRAN Manual, pages 10.7-10.8 for 
details. The Namelist name is $L0AD. 

2 . How to Enter Specific Data Items 

I 

a. Head 

This data set is entered starting in card column 1 and the 
first 42 characters will be picked up. This card must be 
/ the first one in the deck just prior to the Namelist data 
set. 

b. MSELEC 

The total number of rows selected must be < 72. The "from" 
value must be < to the "to" value of each set. The current 
"from" value must be greater than the previous "to". The 
first 0 entry terminates the data set. 

c. INPTRN 

The only valid entries are 5 or 25. 5 - means the data for 

will come via the card reader and 25 - means the data 
MQ 

will be on file = TAPE25. 
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d. 


N0B0D 

Care must be taken to make sure the entry j is between 1 
and NB, where NB is the number of bodies in the model processed 
in UFSSP . j must be a flexible body with non-zero flexible 
degrees of freedom. This quantity must be input. 

e. STIMn (01 < n < 10) 

The current "from” value must be > to the previous "to". 

The "from" value must be < to the "to" value of each set. 

f. HISCSE 

This number must agree with the case number on the UFSSP 
history tape that is mounted. 

C. SAMPLE DATA DECK SETUP 
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V. 0UTPUT 

A. GENERAL 

There is no control on how often data is printed or written on the 
history tape. It is done for each time point for which there is 
a calculation. The history tape is always written. The page is 
restored so that data for an array will not be split between pages. 
Unless otherwise stated, the format of a floating point number is 
+X.XXXXXXE+YY. 

B. PRINTED OUTPUT 

The following types of print exist: 

1. Case header 

2. Error print 

3. Input data print 

4. Summary of L? aX and L™ n - internal loads (maximum and 
minimum values). 

1. Case Header 

This header prints at the top of every page with the following 
format : 

HH — -42 CHARACTER HEADER— HH CASE N0 B0D Y P AGE 

2. Error Print 

See Section VD. for details on types and formats of output. 

Error messages for erroneous input data cards print prior to the 
input data print . 

3. Input Data Print 

This output data set is printed at the beginning of each case 
(on option) and includes all input data quantities and arrays 
as they will be used in the program. All zero values are printed. 
Even if there are input errors which cause the job to be termin- 
ated, this data will be printed. 
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Format : 


a. Floating Point Numbers 

+O.XXXXXXXXE+YY; where trailing zeroes are ignored so that 
the format is variable. 

b'. Fixed Point Numbers 

±XXXX — X; as many positions print out as there are signi- 
ficant digits. Leading zero are ignored. 

c. Single Data Items 

Each individual named data quantity prints one per line as 
follows : 

INPTRN - 5. 

d. Arrays Print as Follows: 

MSELEC - 

FR0M n 2 T0 n 2 FR0M n 3 T0 n 4 FR0M n 5 T0 n g FR0M n ? T0 n g FR0M n g T0 n 1Q 
FR0M 

FR0M n ± _ 3 n i_ 2 FR(2M n i-l n i 
where l<i<50 
STIM = 

FR0M+X . XXXXXXE+YY TO +X . XXXXXXE+YY EVERY n 
FR0M 

FR0M+X . XXXXXXE+YY T0 +X . XXXXXXXE+Y Y EVERY n 
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t t max * T min ^ . 

4. L. ana L. Print 

i 1 

The print format on a page for one L*? aX and L^ ln is as follows: 

MINIMUM L0AD(P0SITION ) = +X.XXXXXXE+YY TMIN=X.XX— 

L0AD VECT0R F0R MINIMUM 


P0S 

VALUE P0S 

VALUE 

P0S 

VALUE 

1 

+X.XXXXXX+YY 2 

+X.XX 

6 

+X.XX 

7 

+X.XX-- 

+X.XX 

12 

+X.XX 


m +X.XX 


MAXIMUM L0AD(P0SITION 

n 

.XXXXXXE+YY TMAX=X . XX_ 

L0AD VECT0R F0R MAXIMUM 



P0S VALUE 

P0S 

VALUE 

1 +X.XXXXXXE+YY 

2 

+X.XX 


The L“ and L^ printout in vector position order 1 < i < m 
where H m i 72. ’ 

C. HISTORY TAPE- 

This tape contains the L^ vector for each time point processed 
by this program. This tape cannot be processed as part of the job 
in which it is created. It is formatted for use by the TRW stand 
alone version of the PL0T8 program. This program will print the 
or generate a plot tape for plotting on a calcomp plotter. See 
the PL0T8 writeup for its use and capabilities. 

1. Format 


I.D. RECORD 


HEAD (1-7) 

- 

"INTERNAL LOADS PROGRAM" 

HEAD (8-14) 

- 

The BCD data from HEAD 

HEAD (15-22) 

- 

Unused 

YEAR. (1-3) 

- 

Unused 

PL0T5 

- 

Unused 

CPLTV (1-48) 


Unused 

TM0N 

- 

Unused 

TM0FF 

- 

Unused 

NWRDS 

- 

Contains m - number of selected rows 

PL0T6 

- 

Unused 

DUMPL8 (1) 

- 

N0CASE 

DUMPL8 (2-30) 

- 

Unused 
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The data which is missing in the above list will be supplied 
via data input to the stand alone PL0T8 program. 



DATA 

RECORD 


TIME 

- 

current time 


L i 

- 

Selector No. 

1 

L 2 

- 

Selector No. 

2 

■4 

- 

Selector No. 

3 

• 

• 

• 


Selector No. 

m 

m 





A data record is m + 1 words long where l<nv<72. These records 
are written in the same way as discussed in Appendices A-E 
of the UFSSP Users Manual, [1] . 

D. ERROR PROCESSING 

1. General 

The program attempts to find all input data card errors during 
input processing that would create havoc with the program logic. 

The remainder of the checks will ensure that the correct loads 
history tape has been mounted and that the physical model it 
represents is compatible with the input card data. The program 
will complete the checking of the input data cards until every 
error has been found. The detection of any error listed below, 
however, will terminate the job so that no calculations will 
be performed. 

2 . List of Errors with Message Printed 
ERROR 

1 A check is made to see if Q/3 from the first card of the 
LTMP loads transformation data is = to NJTS from the 
body dependent header record of the UFSSP history tape. Message 

1 - L0ADS TRANSFORMATION Q IS N0T C0HPATIBLE WITH NJTS 
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2 There is no structures record for flexible body j ' (as found 
via NJPRIM(NOBOD) ) on the loads history tape. This means 
that body j has 0 flexible DOF. Message: 

2 - NOf STRUCTURES DATA RECORD 0N UFSSP HISTORY TAPE F0R B0DY JPRIME 

3 A check is made to make sure the case no on the UFSSP history 
tape is = to HISCSE. Message: 

3 - HISCSE IS N0T EQUAL T 0 UFSSP HIST0RY TAPE CASE N0. L0ADS CASE 

4 The count m of the number of rows selected from the B^ 
matrix via the MSELEC array is > 72 or it is greater than M, 
the number of rows specified in. the first card of the loads 
transformation data. Message: 

4 - T00 MANY R0WS SELECTED FR0M L0ADS TRANSF0RMATI0N MATRIX 

5 Each from-to pair in MSELEC and STIMn arrays are checked to 
make sure that the "FR0M" field is < "T0"‘ field. Message: 

5 - FR0M FIELD IS GREATER THAN THE T0 FIELD. ARRAY ENTRY 

6 Each "FR0M" entry in the MSELEC and STIMn arrays are checked 
to make sure that the current "FROM" entry is > than the 
previous one. Message: 

6 - ENTRIES ARE 0UT 0F 0RDER IN ARRAY 

7 N0B0D is checked to make sure it is in the range 1 to NB. 

NB is obtained from the first record on loads history tape. 
Message : 

7 - N0B0D IS N0T BETWEEN 1 and NB 

8 N0B0D is checked to make sure that NJPRIM(N0B0D) specifies 
a non-zero j'. Message: 

8 - N0B0D SPECIFIES A N0NFLEXIBLE B0DY 

9 INPTRN is checked to be =5 indicating loads transformation 
data is coming via cards or = .25 indicating it is to come 
in via a disk or tape file - TAPE25 . Message: 

9 - INPTRN SPECIFIES AN ERR0NE0US INPUT DEVICE 
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VI. SAMPLE CASE 


This first sample case is presented to illustrate how the required 
input data is to be entered plus the format of the program output. 

A. PHYSICAL DESCRIPTION OF THE CASE 

Figure 1 presents a schematic of the model chosen to demonstrate 
the operation of the program for a simplified configuration. The model 
chosen for simulation is that of a rigid central body with two identical 
flexible beams cantilevered diametrically on the perimeter of the central 
body. Thus, the model consists of three bodies, one rigid and two 
flexible. 

Since the flexible bodies are identical, a single structural 
dynamics model is formulated consisting of a ten-foot long beam with mass 
concentrated at six mass points. Figure 2 depicts the uniform beam model. 
Physical characteristics of the beam are as follows: 

length = 120 in 

/ 2.2 

mass/length = 0.01 lb. sec /in 

El = 6.66 x 10 6 lb. in 2 

yy 7 2 

El = 1.48 x 10' lb. in 

zz 

(The above values for El and El have been chosen to produce a first 
v yy zz 

cantilever bending mode in the x-z plane at a frequency of approximately 
1 Hz and a first cantilever bending mode in the x-y plane at a frequency 
of approximately 1.5 Hz.) As seen in Figure 2, the structural dynamics 
model consists of 9 joints and 7 uniform, massless beam members. The 
total mass of the beam is lumped at joints 2 through 7. Joint 1 is fully 
constrained while joint 8 represents the free end of the cantilever beam 
model. Joint 9 is a fictious joint needed by the structural dynamics 
program to define the three-dimensional geometry of the structure; like 
joint 1, it is fully constrained. 

The model is initially stationary in inertial space. At time 

equal to zero, a triangular force pulse is applied to the center of mass 

1 IP 

of Body 1, directed along the e^ axis (t Q = 0 for F 3 (t) ). In addition, 

at time equal to 2 seconds, a triangular force pulse is applied to the 

1 IP 

Body 1 mass center directed along the e 2 axis (t Q = 2 for F 2 (t) ) . 
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Figure 1. Schematic of Model for Test Case 



The basic objective of this test case is the determination of * 
the interconnection loads between the bodies throughout the time history 
of this simulated thrusting sequence. These loads are computed in two 
ways, thus affording an internal check on the accuracy of the results. 

In the first calculation, the interconnection loads are computed by the 
UFSS Program itself; in the second calculation, the interconnection 
loads are computed by the LOAD Program as being equal to the internal 
loads acting on member 1 at its connection with joint 1 (see Figure 2). 

B. DESCRIPTION OF THE UFSSP Run 

Since the L0AD Program calculates internal loads via information 
obtained from a UFSSP history tape, it is first necessary to generate a 
time history simulation. The following pages contain the Namelist load 
sheets required for this sample case. Since most of the input data is 
the same as that described in the despin test case contained in the UFSSP 
Users Manual, it will not be described in detail here. The only basic 
differences lie in the tabular description of the applied forces. The 
entry ICALCF(l) = 1 indicates that prescribed disturbances are applied to 
Body 1.. ITABN1 (1,3) = 1 and ITABN1(1,2) = 2 further indicate that Body 1 
is acted upon by a prescribed force whose third component is specified 
by the first interpolation table pair and whose second component is specified 
by the second interpolation table pair. 

Figure 3 presents the standard UFSS print out of a typical time 
point (2.5 seconds). The interconnect force and torque values are given 
in units of lb and ft-lb respectively. Their components are along the 

2 o 

je^ body-fixed axes for Body 2 and along the e^ body-fixed axes for Body 3. 

Figures 4 and 5 present plots of the Body 2 interconnection 

forces and torques respectively. Note that F^ 1 and T? 1 are due to the 
IP 3 2 

applied force F 3 while the remaining interconnection forces and torques 

arise^from the applied force F^. The values for F^ 1 are very small 

10 ) as they should be from symmetry considerations. The values for 

2 1 

■^1 are small and arise purely through the flexible deformations of the 
beam. 
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C. DESCRIPTION OF THE LOAD RUN 


For this- test case the member forces at the start of member 1 
are calculated via the appropriate LTM and Equation (1) for Body 2. 

Since Joint 1 has zero mass, these member forces should correspond to 
the interconnection forces and torques as calculated in the UFSS run 
described in B. above. This present section will discuss details of the 
L0AD run only; for details regarding definition of the internal loads 
and the LTM itself, the user is referred to Appendix A. 

The Namelist input data set for this run is as follows: 

THRUSTER L0AD TEST CASE 

$L0AD 


N0CASE 

= 

002, 

N0B0D 

= 

2, 

INPTRN 

= 

25, 

MSELEC 

(1,1) = 1,6, 

HISCSE 

= • 

400, 

STIM01 (1) 

- 0,10,1, 

INPRNT 

$END 


1, 


The N0CASE entry merely serves to identify the run. The "N0B0D = 2" 
entry indicates that internal loads are to be computed for Body 2, while 
the INPTRN entry specifies that the loads transformation matrix is to be 
input via a disk file. The MSELEC entry indicates that rows one through 
six of the LTM are to read into core for calculation, while the "HISCSE = 
400" entry is set to match the "ICASE = 400" entry of the UFSS run used 
to supply input data; this correspondence assures that the correct UFSS 
history tape is read. 

The entry "STIMOl(l) = 0,10,1" indicates that loads are to be 
computed for every time point on the UFSS history tape (recall that the 
UFSS run started at time equal to zero and ended at time equal to ten 
seconds). Finally, the "INPRNT = 1" entry specifies that the input card 
data set should be printed. 
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Figure 6 presents the standard print out of the input card 
data set. The print out is self-explanatory. 


Figure 7 presents the standard output print of the LOAD Program. 
The correspondence between the load positions (rows of the LTM) and the 
interconnection forces calculated by the UFSS Program are as follows: 

(See Appendix A) 


Load Position 
1 

2 

3 

4 

5 

6 


UFSS Quantity 




F 

T 

T 

T 


21 

3 

21 

1 

21 

2 

21 

3 


Thus, for example. Figure 7 reveals that F^ attains its maximum value of 
5.97228 lb at four seconds into the run while attaining its minimum value 
of -4.591991 lb at 0.95 seconds into the run. All of the maxima and 
minima presented in Figure 7 can readily be verified by inspection of the 
plots for the member loads presented in Figure 8 (these plots are obtained 
by using the TRW stand-alone program AC041, with the LOAD history tape as 
input data). It should be noted here that the LTM was produced by the 
TRW program "LTMP" which — like the TRW program, "SMAP", used to generate 
the structures data for the test case — accepts only the units of lb, 
inch and second. Thus, all member loads corresponding to moments are 
output in units of lb-in by the LOAD Program. Finally, the following 
comparison of load values, as computed by the UFSS Program and the LOAD 
Program, are presented for a typical time point (2.5 seconds): (Units 

for moments are lb. in) 
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UFSSP 


LOAD 



1.34644x10 ^ 

1.346442x10 


3.87862 

3.878621 


2.76538 

2.765376 


1 . 9054x1 0 _1 

0 


-2.4074xl0 2 

-2.407526xl0 2 


3. 3765xl0 2 

3.376713xlQ 2 


It should be noted that the value of as predicted by the LOAD 
Program is identically zero. This computation is inherent to the 
geometry of the beam; since all mass points lie on the e£ axis, there 
is no moment arm to allow either applied or d'Alembert forces to produce 
a moment about the e^ ax i s . 

It should be noted that care must be exercised when computing internal 
loads for. a rapidly spinning flexible body. For the case of the model shown 
in Figure 1, computation of internal loads for a steady spin about the e^ 
axis will yield incorrect results for the member moments. These moments 
contain contributions due to the product of the centrifugal d'Alembert 
force and the instantaneous flexible deformations in the e^ and e^ directions 
at each joint. Since the load transformation matrix used in Equation (1) 
is a constant matrix, there is no way that these moments can be included. 
Calculations for the member forces are correct. 

The calculations of the interconnection forces and moments as performed 
within the UFSSP simulation run are correct for any motion of the model. 
Therefore, it is always advisable that these UFSSP values be compared to the 
root member forces computed by the L0AD Program in order that such spin 
effects may be identified and, if necessary, added to the computed member 
forces. 
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Figure 6. Standard Print-out o£ LOAD Input Data . 








Figure 7. LOAD Standard Output Print 
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PART II. OPERATION MANUAL 


VII. SOFTWARE /HARDWARE REQUIREMENTS 
CPC 6500 

This program runs under the TRW / TS S system and requires 41,000 words of 
• core and 3 files. These files may be on disk or tape. 

UHIVAC 1108 

This program runs under the exec 8 system and requires 40,000 words of 
core plus 3 files. These files can be on drum or tape. 
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VIII. OVERALL DECK SETUP 


CPC 6500 

The following diagram assumes that the program, in relocatible form, is 

stored on disk under file name L0ADPG and the matrix is to be 

MO 

input from cards. 
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UN I VAC 1108 


(To be supplied by MSFC programming staff 
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IX. tape/disk/drum usage 


Program 

Symbol 

CDC 

File 

No. 

UNIVAC 

File 

No. 

Tape 

Usage 

Mode 

INPTRN 

25 

25 

bA. Matrix Card File 
MQ 

BCD 

10 

10 

10 

UFSSP History Tape 

BINARY 

8 

8 

8 

Loads History Tape for Plotting 

BINARY 

5 

5 

5 

Input Card Set 

BCD 

6 

6 

6 

Print Data Set 

BCD 
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XI 


BLOCK DIAGRAM OF LOAD PROGRAM 



END 
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Appendix A. Load Transformation Matrix Description 


This appendix further defines and describes the flexible body 
structural dynamics model and the internal member loads within a structure 
modeled as a system of lumped-mass joints which are interconnected by 
weightless beam elements. The following discussion is mainly drawn from 
the user manual for the TRW developed LOAD TRANSFORMATION MATRIX PROGRAM 
(LTMP) — the program used to generate the LTM for the sample case 
contained in Section VI of this document. For purposes of analysis, the 
structures admissible to the LTMP are grouped into the following five 
types: 

1) Three-dimensional trusses, 

2) Three-dimensional frames whose member cross sections 
have polar symmetry, 

3) Two-dimensional frames (motion in XY plane of structure), 

4) Two-dimensional grids (motion normal to XY plane of structure), 

5) Three-dimensional frames whose member cross sections need 
not have polar symmetry. 

For use with the UFSS and LOAD Programs, we will consider only structure 
type 5. 

Typical of most structural dynamics programs, the following 
quantities describe the structure: 

1) Coordinates of the joints, 

2) Properties of the materials of which the structure is 
composed, 

3) Geometric properties of the member elements, 

4) Location of restraints. 

The elastic properties for the elements are defined by the material 
properties. The material properties consist of Young’s modulus of 
elasticity, Poisson's ratio, the shear modulus, and the coefficient of 
thermal expansion in each of two orthogonal directions (the local x and 
y axes of the element). The material density is also entered. 
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The geometric properties of a member element consist of the joint 
numbers at the start and end of the member, the cross sectional area, 
the effective shear areas in the local y and z axes, the area moment of 
inertia about the local y and z axes, and the section torsional stiffness. 

Two types of reference systems are used in the LTMP : right-handed, 

orthogonal, Cartesian coordinate systems for structure geometry, displace- 
ments, and forces; and a numbering convention to define the degrees of 
freedom in the structure. 

Coordinate Reference Systems 

Right-handed, orthogonal, Cartesian coordinate systems are used 
to describe the geometry of the structure, the joint displacements, 
and the member displacements and forces. The description of the 
structure requires two different coordinate systems. One is the 
overall system, termed the global coordinate system. The other is 
a coordinate system which is fixed to each element in the structure; 
these latter are called local coordinate systems. 

Global Coordinate System 

The global coordinate system is an arbitrarily selected coordinate 
system, usually positioned so the directions of the axes coincide 
with some of the major dimensions of the structure. Upper case 
letters are used when referring to the global coordinate system. 

The positive directions are denoted by X, Y, Z. (These directions 
coincide with the body-fixed axes _e^, e^ and e^ used in the UFSS 
and LOAD Programs.) The global coordinate system is used to 
describe joint coordinates, joint deflections, and externally applied 
forces. The six possible displacement components (degrees of freedom), 

fi jX» 5 jy» ®j X » y ’ 9 jZ’ f ° r the 3 th j° int are defined in the 

global coordinate system as shown in Figure A.l. (These displace- 
ment components are identical to the joint modal deformations <j>^ N , 

^2 * ^3 * 9 1 ’ 9 2 anc ^ 9 3 ^ or N of Body j as used in the 

UFSS and LOAD Programs.) 
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Reference System for Degrees of Freedom 

A number is assigned to each degree of freedom in the structure, 
whether it is rigidly constrained or free to move. The degrees of 
freedom are numbered sequentially, starting at joint 1 and 
continuing to the last joint. Thus, for example, the degrees of 
freedom at joint 5 of a structure would be numbered as shown below. 


V 


25 



26 

6 jz 


27 

V 

• 

28 

V 

• 

29 

6 JZ 

• 

30 


Member Element 

The idealized member elements used to represent a weightless beam 
or rod between two points in the structure are subject to the 
following assumptions: 

1) All members are represented by a straight line which is 
attached to two joints of the structure. The straight 
line corresponds to the centroidal axis of the member. 
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2) All members have uniform cross sectional properties between 
the two end joints with the section properties defined about 
the principal axes of the cross section. 

3) The local y and z axes are the principal axes of the cross 
section. 

4) The shear center of the cross section coincides with the 
centroid of the section. 

5) Axial forces act along the centroidal axis of the member. 

6) The member's elemental stiffness matrix is formed using 
small deflection theory. 

Since the member elements are interconnected at discrete points, 
the elemental stiffness matrix of the member represents the exact 
stiffness relationship between the two end joints. 

Local Coordinate System 

The position and orientation of a given member element in space is 
specified by the coordinates of the end joints and a linear trans- 
formation matrix relating the local member coordinate system to the 
global coordinate system. The local coordinate system is fixed to 
the member and is used to relate the stiffness properties, deflections, 
and forces of the member back to the global coordinate system. Lower 
case letters are used when referring to the local coordinate systems. 
The local x axis always extends in a positive direction from the 
start of the member (first joint specified) to the end of the member 
(second joint specified). The orientation of the y and z axes is 
determined as follows: 

Type 5 . A reference joint (specified by input data) is required 
to define the orientation of the y and z axes for the Type 5 
structure. Let R be a vector from the first joint of the member 
to the reference joint (see Figure A. 2). The local z axis is 
perpendicular to the plane formed by the member's local x axis 
and the vector R. It is positive in the direction of the vector 
cross-product of the local x axis and R respectively. The local 
y axis is in the plane which is formed and in a direction which 
will form a right-handed, orthogonal, Cartesian coordinate system 
with the local x and z axes. Note that the reference joint must 
not be located on the local x axis. The reference joints must be 
an existing joint in the input data. If none of the actual joints 
in the structure can be used as a reference joint, an arbitrary 
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Type 5 (Cont'd) 

reference joint with all degrees of freedom rigidly constrained 
can be specified without increasing the size of the stiffness 
matrix. 

Member Displacements and Forces 

Each member has twelve possible displacements defined in the local 
coordinate system. These displacements are obtained by transforming 
the joint displacements at the ends of the member from the global 
coordinate system to the local coordinate system of the member. A 
member displacement is positive in the positive direction of the 
local coordinate system. 

The member forces are also defined in the local member coordinate 
system. Member forces are positive if they act in the positive 
direction of the member coordinate system. The twelve possible 
forces associated with a member of the structure are shown in 
Figure A. 3. All forces shown are positive. 



Figure A. 3. Member Forces 
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Load Transformation Matrix Definition 


As defined previously, each row of the LTM corresponds to a given 
member force. Thus, there can be a total of 12M member forces 
(internal loads) computed, where M is the total number of members 
within the structural dynamics model. 

Let us now consider the structural model used in the test case of 
Section VI. The LTM generated for this structure is presented in 
Figure A. 4. The first line of the card file indicates that the 
LTM contains 12 rows and 27 columns. The number of columns is 
completely determined by the model as being equal to three times 
the number of joints. However, the number of rows is arbitrary, 
being equal to the number of member forces (internal loads) selected 
by the user. For the present case, the load coefficients for the 
six member forces at the start of member 1 (i.e., at its connection 
with joint 1) comprise the elements of the first six rows of the 
LTM. The load coefficients for the six member forces at the start 
of member 5 (i.e., at its connection with joint 5) comprise the 
elements of the last six rows of the LTM. 

It should be noted that columns 1, 2 and 3 (corresponding to 
external forces applied to joint 1) and columns 25, 26 and 27 
(corresponding to external forces applied to joint 9) are identically 
zero due to the fact that these two joints are fully constrained. 

In addition, rows 4 and 10 (corresponding to axial moments within 
members 1 and 5) are identically zero due to the fact that all 
joint masses lie on the global X axis which corresponds to the local 
x axis for each member. 

Through the MSELEC entry in the sample case, only the member forces 
in member 1 were calculated in the sample case run. 
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Figure A. 4. Load Transformation Matrix for Sample Case (continued) 
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